One particle spectral weight of the three dimensional single band Hubbard model 
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Dynamic properties of the three-dimensional single-band Hubbard model are studied using Quantum 
Monte Carlo combined with the maximum entropy technique. At half-filling, there is a clear gap in 
the density of states and well-defined quasiparticle peaks at the top (bottom) of the lower (upper) 
Hubbard band. We find an antiferromagnetically induced weight above the naive Fermi momentum. 
Upon hole doping, the chemical potential /j, moves to the top of the lower band where a robust peak 
is observed. Results are compared with spin-density-wave (SDW) mean-field and self consistent 
Born approximation results, and also with the infinite dimensional (D — oo) Hubbard model, and 
experimental photoemission (PES) for three dimensional transition-metal oxides. 
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I. INTRODUCTION 

The single band two dimensional Hubbard Hamilto- 
nian Q has recently received considerable attention due 
to possible connections with high temperature super- 
conductors. Indeed, evidence is accumulating that this 
Hamiltonian may describe, at least qualitatively, some of 
the normal state properties of the cuprates. Exact Di- 
agonalization (ED) and Quantum Monte Carlo (QMC) 
have been used to model static properties like the be- 
havior of spin correlations and magnetic susceptibility 
both at half-filling and with doping. Q Comparisons of 
dynamic quantities like the spectral weight and density 
of states with angle-resolved photoemission results 
have also proven quite successful. Significantly, while an- 
alytic calculations have pointed towards various low tem- 
perature superconducting instabilities, such indications 
have been absent in numerical work. |2j 

Historically, however, the Hubbard model was first pro- 
posed to model magnetism and metal-insulator transi- 
tions in 3D transition metals and their oxides, |l]] rather 
than superconductivity. Now that the technology of nu- 
merical work has developed, it is useful to reconsider 
some of these original problems. A discussion of possible 
links between the 3D Hubbard model and photoemission 
results for YTi0 3 , SrV0 3 and others §-0 has already 
recently occurred. In such perovskite Ti""~and V 4+ ox- 
ides, which are both in a 3d 1 configuration, the hopping 
amplitude t between transition-metal ions can be varied 
by modifying the d — d neighboring overlaps through a 
tetragonal distortion. Thus, the strength of the electron 
correlation U/t can be varied by changing the composi- 
tion. In fact, a metal-insulator transition has been re- 
ported in the series SrV0 3 - CaV0 3 - LaTi0 3 - YTi0 3 . 
On the metallic side, a quasiparticle band is experimen- 
tally observed near the Fermi energy Ep , as well as a 
high energy satellite associated to the lower Hubbard 
band (LHB). pj|l6|| Spectral weight is transferred from 
the quasiparticle to the LHB as U/t is increased at half- 
filling. 

In this paper, we report the first use of Quantum 



Monte Carlo, combined with analytic continuation tech- 
niques, to evaluate the spectral function and density of 
states for the 3D Hubbard Hamiltonian. The motivation 
is twofold. First, we want to compare general proper- 
ties of the 3D Hubbard Hamiltonian with the extensive 
studies already reported in the literature |jll|,^2|,^6|j|] for 
the 2D and infmite-D cases. Of particular importance is 
the presence of quasiparticles near the half-filling regime, 
as well as the evolution of spectral weight with doping. 
Many of the high-Tc cuprates contain Cu02 planes that 
are at least weakly coupled with each other, and thus the 
study of the 3D system may help in understanding part 
of the details of the cuprates. More generally, the Hub- 
bard Hamiltonian is likely to continue being one of the 
models used to capture the physics of strongly correlated 
electrons, so we believe it is important to document its 
properties in as many environments as possible for po- 
tential future comparisons against experiments. 

Secondly, we discuss a particular illustration of such 
contact between Hubbard Hamiltonian physics and ex- 
periments on 3D transition metal oxides. In addition 
to the studies of half-filled systems with varying corre- 
lation energy mentioned above, experiments where the 
band filling is tuned by changing the chemical composi- 
tion have also been reported. |T^,|lC],|l4| One compound 
that has been carefully investigated in this context is 
Yi_ x Ca x Ti0 3 . At x = the system is an antiferromag- 
netic insulator. As x increases, a metal-insulator tran- 
sition is observed in PES studies. The lower and up- 
per Hubbard bands (LHB and UHB) are easily identified 
even with x close to 1, which would naively correspond 
to small electronic density in the single band Hubbard 
model, i.e. a regime where U/t is mostly irrelevant. In 
the experiments, a very small amount of spectral weight 
is transferred to the Fermi energy, filling the gap observed 
at half- filling (i.e. generating a "pseudogap" ) . 

Analysis of the PES results of these compounds using 
the paramagnetic solution of the Hubbard Hamiltonian in 
infmite-D p5| , a limit where dynamic mean field theory 
becomes exact (see section II), has resulted in qualita- 
tive agreement [0,|l7],|l6| with the experimental results. 
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At and close to half-filling there is an antiferromagnetic 
(AF) solution which becomes unstable against a param- 
agnetic (PM) solution at a critical concentration of holes. 
In the PM case, weight appears in the original Hubbard 
gap as reported experimentally. However, this analysis of 
the spectral weight in terms of the infinite-D Hamiltonian 
is in contradiction with results for the density of states 
reported in the 2D Hubbard model Q where it is found 
that upon hole (electron) doping away from half-filling 
the chemical potential \i moves to the top (bottom) of 
the valence (conduction) band. The results at (n) = 1 
in 2D already show the presence of a robust quasiparti- 
cle peak which is absent in the insulating PM solution 
of the D = oo model. That is, in the 2D system the 
large peak in the density of states observed away from 
half-filling seems to evolve from a robust peak already 
present at half-filling. On the other hand, at D = oo 
a feature resembling a "Kondo-resonance" is generated 
upon doping if the paramagnetic solution is used. This 
peak in the density of states does not have an analog 
at half-filling unless frustration is included. Q Studies 
in 3D may help in the resolution of this apparent non- 
continuity of the physics of the Hubbard model when the 
dimension changes from 2 to oo. The proper way to carry 
out a comparison between D = 3 and oo features is to 
base the analysis on ground state properties. With this 
restriction, i.e. using the AF solution at D — oo and close 
to half-filling, rather than the PM solution, we found that 
the D = 3 and oo results are in good agreement. 

In this paper we will consider which of these situations 
the 3D Hubbard Hamiltonian better corresponds to, and 
therefore whether the single band Hubbard Hamiltonian 
provides an adequate description of the density of states 
of 3D transition- metal oxides. 



II. MODEL AND METHODS 

The single band Hubbard Hamiltonian is 
H = -t^{c\ a c- }a + h.c.) - fiJ2n ia 

(ij> i<J 

+[/^> iT -l/2)(n u -l/2), (1) 

i 

where the notation is standard. Here (ij) represents 
nearest-neighbor links on a 3D cubic lattice. The chemi- 
cal potential u controls the doping. For fj, = the system 
is at half filling ((n) = 1) due to particle- hole symmetry. 
t = 1 will set our energy scale. 

We will study the 3D Hubbard Hamiltonian using a fi- 
nite temperature, grand canonical Quantum Monte Carlo 
(QMC) method which is stabilized at low temper- 
atures by the use of orthogonalization techniques [pL9|| . 
The algorithm is based on a functional-integral represen- 
tation of the partition function obtained by discretizing 
the "imaginary-time" interval [0, 0] where is the inverse 



temperature. The Hubbard interaction is decoupled by 
a two-valued Hubbard-Stratonovich transformation |2Q] 
yielding a bilinear time-dependent fermionic action. The 
fcrmionic degrees of freedom can hence be integrated out 
analytically, and the partition function (as well as ob- 
servables) can be written as a sum over the auxiliary 
fields with a weight proportional to the product of two 
determinants, one for each spin species. At half-filling 
((n) = 1), it can be shown by particle-hole transforma- 
tion of one spin species (q| — > (— l)*Cj|) that the two 
determinants differ only by a positive factor, hence their 
product is positive definite. At general fillings, however, 
the product can become negative, and this "minus-sign 
problem" restricts the application of QMC to relatively 
high temperature (of order 1/30 of the bandwidth) off 
half-filling. 

The QMC algorithm provides a variety of static and 
dynamic observables. One equal time quantity in which 
we are interested is the magnetic (spin-spin) correlation 
function, 

j 

Here m- } = ^ CT an^ a is the local spin operator, and N is 
the total number of lattice sites. Static correlations have 
also been investigated in earlier studies of the 3D Hub- 
bard model (2l].|z| where the antiferromagnetic phase di- 
agram at half filling was explored. 

To obtain dynamical quantities in real time or fre- 
quency, the QMC results in imaginary time have to be 
analytically continued to the real time axis. Since we 
are mostly interested in the one-particle spectrum we 
measure the one-particle Green function G(p,r). The 
imaginary part of G(p,uj) (in real frequency) defines the 
spectral weight function at momentum p, A(p, lu), which 
is related to G(p, r) by: 

/oo „-ra 
duAfau)——^. (3) 
-oo i- -\- e 

A(p,u>) can in principle be calculated by inverting Eq. (||) , 
but the exponential behavior of the kernel at large values 
of \u>\ makes this inversion difficult numerically. G(p,r) 
is quite insensitive to details of A(p, lu) in particular at 
large frequencies. Since G(p, r) is known only on a finite 
grid in the interval [0, 0] and there only within the statis- 
tical errors given by the QMC-sampling, solving Eq.(||) 
for A(p, lo) is an ill-posed problem. A large number of 
solutions exists, and the problem is to find criteria to se- 
lect out the correct one. This can be done by employing 
the Maximum Entropy (ME) method (23). Basically, the 
ME finds the "most likely" solution A(p, ui) which is con- 
sistent with the data and all information that is known 
about the solution (like positivity, normalization, etc.). 
ME avoids "overfitting" to the data by a "smoothing" 
technique that tries to assimilate the resulting A(p, lu) to 
a flat default model. In the absence of any data (G(p, r)) 
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ME would converge to the default model which is chosen 
to be a constant within some large frequency interval. 
There is no adjustable parameter in the ME application. 
One needs accurate data for G(p,r) with a statistical 
error of o(10 -4 ) to get reliable results for A(p,w). 

In principle, one can calculate analytically the first and 
second (and higher) moments of the spectral weight and 
include this information in the ME procedure. However, 
we chose to calculate the moments afterwards from the 
resulting function A(p, u>), and to compare them with the 
analytically known results as a further test. The agree- 
ment was in all cases within 10%. Still, the ME meth- 
ods provides only a rough estimate of the true spectral 
weight functions. Band gaps and the positions of signif- 
icant peaks are usually well captured but fine structure 
which needs a high frequency resolution is hard to detect 
within the ME approach. 

Integrating A(p, lu) over the momenta gives the one- 
particle density of states (DOS) N(u). However, tech- 
nically, it is preferrable to integrate first G(p, r), which 
reduces the statistical errors, and then perform the ana- 
lytic continuation. 

The DOS will be compared to results from the dynam- 
ical mean- field theory of infinite dimensions, D = oo, 
JL5| . In this limit, with the proper scaling of the hop- 
ping element (t = t*/\fZ, with Z being the coordina- 
tion number) the one-particle self energy becomes local 
or, equivalently, momentum independent and the lattice 
problem is mapped onto a single site problem. The con- 
stant t* is set to t* — \/6 to obtain the same energy scale, 
t = 1, when compared to the 3D case. J24j In contrast to 
conventional mean- field theories, the self energy remains 
frequency dependent, preserving important physics. Spa- 
tial fluctuations are neglected, an approximation which 
becomes exact in the limit Z — > oo (Z — 2D for the 
simple cubic lattice). Even in D = oo, the remaining 
local interacting problem cannot be solved analytically 
but will also be treated by a finite temperature QMC 
p5[ supplemented by a self-consistency iteration p7|,^6[ . 
The advantage is that the system can be investigated in 
the thermodynamic limit with a modest amount of com- 
puter time. Due to its local character the D = oo ap- 
proach cannot provide information on momentum depen- 
dent spectral functions. However, recently a /c-resolved 
spectral function has been studied (2(| in D = oo. 

Among other things, the D — oo limit has been used 
recently to study the AF phase diagram in the Hubbard 
model |l2| . The agreement of the Neel temperature with 
3D results is good. Q In D = oo it is further possible 
to suppress AF long range order artificially by restrict- 
ing the calculation to the (at low temperatures unsta- 
ble) paramagnetic solution at half-filling. In this way, 
one may simulate frustration due to the lattice structure 
or orbital degeneracy, although in the absence of cal- 
culations for hypercubic lattices with nearest and next- 
nearest uniform hopping amplitudes it is still a conjecture 
how close this approach is to including these effects fully. 



III. HALF-FILLING 

A. Quantum Monte Carlo 

We first study the single particle spectral weight 
A(p,uj) at relatively strong coupling, U = 8, and half- 
filling ((n) = 1) at a low temperature of T — 1/10. 




00 

FIG. 1. Single particle spectral function A(p,u) versus u> 
for several momenta. The results correspond to a lattice with 
4 3 sites, U = 8, P = 10, half-filling and using i = 1 as energy 
scale. Momenta, p, are in units of n/2. Bars indicate the 
position of the quasiparticle peak. 

A gap is clearly present in the spectrum (Fig. 1) which 
is compatible with the expectation that the half-filled 
Hubbard model on a bipartite lattice is an antiferromag- 
netic insulator for all nonzero values of the coupling U/t. 
The spectral weight has four distinct features (two in the 
LHB and two identical ones in the UHB, as expected from 
particle-hole symmetry). In the UHB there is weight at 
a high energy, roughly in the interval between ~ 5t and 
8t. This broad feature likely corresponds to the incoher- 
ent part of the spectral function found in previous sim- 
ulations for the 2D Hubbard and t — J models. M The 
dominant scale of this incoherent weight is t, and since it 
is located far from the top of the valence band its pres- 
ence is not important for the low temperature properties 
of the system. 

Much more interesting is the sharper peak found close 
to the gap in the spectrum. This band dispersion starts 
at a binding energy of approximately lu — —At at mo- 
menta (0, 0, 0) and moves up in energy obtaining its max- 
imum value at lu ss — 2t at momenta (0, 7r/2,7r) and 
(it/ 2, 7r/2, 7t/2) in Fig. 1. The width of the peak dimin- 
ishes as the top of the valence band is reached. Similar 
structure was discussed before in studies of 2D systems, 
which had a somewhat higher resolution, as a "quasipar- 
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tide" band corresponding to a hole moving coherently in 
an antiferromagnetic background. J^,^7j This quasiparti- 
cle should be visualized as a hole distorting the AF order 
parameter in its vicinity. In this respect it is like a spin- 
polaron or spin-bag, [E8) although "string states" likely 
influence its dispersion and shape. 0] The quasiparticle 
(hole plus spin distortion) movement is regulated by the 
exchange </, rather than t. 

Using the center of the quasiparticle peaks of Fig. 1 as 
an indication of the actual quasiparticle pole position, we 
obtain a bandwidth W of about 2 to 3£ or, equivalcntly, 
4 to 6 J using J = 4t 2 /U for U = 8. However, due to 
the low resolution of the ME procedure, reflected in part 
in the large width of the peaks of Fig. 1, it is difficult 
to show more convincingly within QMC/ME that the 
quasiparticle bandwidth is indeed dominated by J. 

Note that moving from (0,0,0) to (ir,ir,ir) along the 
main diagonal of the Brillouin Zone (BZ), the PES part 
of the spectrum (i.e. the weight at uj < 0) loses intensity. 
There is a clear transfer of weight from PES at small 
|p| to IPES at large |p|, as observed in 2D simulations. 
p9| In addition, note that there is PES weight above the 
(naive) Fermi momentum of this half-filled system. For 
example, at p = (0,7r,7r), spectral weight at u> < can 
be clearly observed. Similarly, at p = (0, 0, ir) weight 
in the IPES region is found. This effective doubling of 
the size of the unit cell in all three directions is a con- 
sequence of the presence of AF long range order. The 
hole energy at p = (0, 0, 0) and (n, ir, tt) becomes de- 
generate in the bulk limit and the quasiparticle band, 
for example along the main diagonal of the BZ, has a 
reflection symmetry with respect to (7r/2, ir/2, 7r/2), as 
observed in our results (Fig. 1). However, note that the 
actual intensity of the AF-induced PES weight close to 
(tt, 7r,7r) is a function of the coupling. As U/t — > 0, 
the intensity of the AF-induced region is also reduced 
to zero. The presence of this AF-generated feature has 
recently received attention in the context of the 2D high- 
Tc cuprates. [^8 30 While its presence in PES experi- 
ments at optimal doping is still under discussion, these 
features clearly appear in PES experimental studies of 
half-filled insulators, like Sr 2 Cu0 2 Cl 2 . Q Thus, while 
this behavior has been primarily discussed in the context 
of 2D systems angle-resolved PES (ARPES) studies of 
3D insulators like LaTiC>3 might also show such features. 

In Fig. 2 we show the density of states (DOS) N(ui) of 
the 4 3 lattice. The two features described before, namely 
quasiparticle and incoherent background, in both PES 
and IPES are clearly visible. Also shown in Fig. 2 is the 
temperature effect on N(u>) which is weak for the given 
temperatures ((3 = 10 and 4). The basic features are still 
retained, only the gap is slightly reduced and the quasi- 
particle peak less pronounced at the higher temperature. 

Results corresponding to a larger coupling, U =12, 
at f3 = 4 are shown in Fig. 3. The gap increases and 
the quasiparticle band becomes sharper as U /t grows, as 
expected if its bandwidth is regulated by J. A(p, uj) at 
U = 12, (3 = 4 are similar to the results shown in Fig. 1. 
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FIG. 2. Density of states N(ui) of the 3D Hubbard model 
on a 4 3 lattice, U = 8, (3 = 10 (solid line) and (3 = 4 (dashed 
line) and (n) = 1. We do not enforce the p{u>) = p(— uj) sym- 
metry which occurs at half-filling due to particle-hole sym- 
metry. However, deviations from this constraint are small, a 
further check on our numerics. 
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FIG. 3. Density of states N(u>) of the 3D Hubbard model 
on a 4 3 lattice at U — 8 (dashed line) and U = 12 (solid line), 
f3 = 4, and (n) = 1. 

A characteristic double peak like that seen in Figs. 2,3 
has been observed in the X-ray absorption spectrum of 
LaFeC>3 (S^j, which is a strongly correlated, wide gapped 
antiferromagnetic insulator. The peaks appear at ener- 
gies of about 2.2 eV and 3.8 eV. When Fe is substituted 
by Ni this structure vanishes and the system becomes a 
paramagnetic metal at a Ni concentration of about 80%. 
If we choose for comparison with the calculated DOS 
(Fig. 3) a hopping amplitude of t = 0.5 eV, giving a 
reasonable d-bandwidth of W = 6 eV, the positions of 
the quasiparticle peak and the maximum of the incoher- 
ent band for U = 12t are at about oji « At — 2 eV 



4 



and u>2 ~ 8 At = 4.2 eV, respectively. The agreement 
with the experimental values is fairly good considering 
the crude simplifications of the Hubbard model such as 
neglecting orbital degeneracy and charge transfer effects. 
Even the estimated charge gap of Fig. 3, defined by the 
onset of spectral weight relative to the Fermi energy, 
^charge ~ 3t = 1.5 eV is not too far from the experi- 
mental value of ~ 1.1 eV. The ratio r — uji/loi decreases 
with U since for U ^> t both energies are expected to 
converge to U/2. While for U = I2t, r w 2.1 is compara- 
ble to the experimental value (~ 1.7), it is too large for 
U = 8 (r » 2.9) showing that under the assumption of a 
single band Hubbard model description for LaFeC>3, the 
effective on-site interaction is at least of the size of the 
d-bandwidth. 

Another feature which has been attributed to antifer- 
romagnetic ordering was found in a high resolution PES 
study of V 2 3 . [HQ In the AF insulator at T = 100K 
the spectrum shows a shoulder at cjj = —0.8 eV which 
is absent in the paramagnetic metal at T — 200K . This 
shoulder might be a reminiscent of the quasiparticle peak. 
The maximum of the lower Hubbard band is at about 
LU2 ~ —1.3 eV, giving a ratio 0*2 /wi ~ 1.6 similar to 
that observed in LaFeC>3. The on-site interaction was 
estimated to be about 1.5 times the bandwidth. [B3| 
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FIG. 4. Density of states N(u) of the 3D Hubbard model 
on a 4 3 lattice at U = 12, (3 = 4 and (n) = 1 (solid line) 
compared with the density of states at the same parameters 
for D — 00 (dashed line). 

It is interesting to compare the results obtained in 
our simulations with those found in the D = 00 limit 
of the Hubbard model. At half-filling for arbitrary cou- 
pling strength, the D — 00 model has an AF insulating 
ground state. Its DOS is shown in Fig. 4, using the same 
coupling and temperature as in the 3D simulation. N(uj) 
for both cases are similar, and they are also similar to 
results found before in 2D, suggesting that the physics of 
holes in an antiferromagnetic system is qualitatively the 
same irrespective of whether a 2D, 3D or 00D lattice is 



used, at least within the accuracy of present QMC/ME 
simulations. 



B. SDW mean-field and Born approximation 

Since the data shown in the previous subsection corre- 
spond to holes in a system with AF long-range order, it 
is natural to compare our results against those found in 
mean-field approximations to the half-filled 3D Hubbard 
model that incorporate magnetic order in the ground 
state. The "spin-density-wave" mean-field approxima- 
tion has been extensively used in the context of the 2D 
Hubbard model, fl35| and here we will apply it to our 
3D problem. For a lattice of N sites, the self-consistent 
equation for the gap A is 
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2N 



Si 



(4) 



where Ey, 



j p — y r P + A 2 i s the quasiparticle energy, and 

e p = — 2t(cos p x + cos p y + cos p z ) is the bare electron dis- 
persion. The resulting quasiparticle dispersion is shown 
in Fig. 5 compared against the results of the QMC/ME 
simulation. The overall agreement is good if the cou- 
pling U in the gap equation Eq.(4) is tuned to a value 
U ~ 5.6. It is reasonable that a reduced U should be 
required for such a fit, since the SDW MF gap is usu- 
ally larger than the more accurate QMC result. Similar 
renormalizations of U in comparing QMC and approx- 
imate analytic work has been discussed in the context 
of fitting the magnetic response, [|(| and has also been 
explicitly calculated. |37| Fig. 5 shows many of the fea- 
tures observed in the numerical simulation, namely a hole 
dispersion which is maximized at (ir/2,ir/2,ir/2) for the 
momenta shown there, an overall bandwidth smaller than 
the noninteracting one, and the presence of AF-induced 
features in the dispersion above the naive Fermi momen- 
tum. 

Thus the SDW MF approach qualitatively captures 
the correct hole quasiparticle bandwidth J at half-filling. 
However, a spurious degeneracy appears in the hole 
dispersion in this approximation. Momenta satisfying 
cos + cos p v + cos p z = have the same energy. This is 
not induced by symmetry arguments and is an artifact of 
the SDW approach. In addition, A(p, uj) in the SDW ap- 
proximation only has one peak in the PES region for each 
value of the momentum, missing entirely the incoherent 
part. 

While it may not be necessary to fix this problem in 
this case, it is important in general to be able to go be- 
yond SDW MF. To do this, the self-consistent Born ap- 
proximation HU (SCBA) for one hole in the 3D t — J 
model, which corresponds to the strong coupling limit of 
the Hubbard model, can be used. This technique repro- 
duces accurately exact diagonalization results in the 2D 
case. [§8| 
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FIG. 5. Quasiparticle dispersion in the SDW MF approxi- 
mation at half-filling and zero temperature. Results are at 
U/t — 8 (dot-dashed line) and U = 5.57i (dashed line). 
QMC/ME results for the 3D Hubbard model on a 4 3 lattice 
at the same coupling, density and temperature are also shown 
(open circles). The area in the circles is proportional to the 
peak intensity. The error bars are the width of the peak. For 
some momenta the intensity of the PES or IPES data is so 
low that no peak position is reported. 
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FIG. 6. Quasiparticle dispersion in the SDW MF approxi- 
mation at U = 12. 8t (dashed line) compared against similar 
results obtained with the SCBA approximation for one hole 
in the t — J model at J/t = 0.3125 (solid line). The energy 
scale of the SCBA dispersion is shifted such that the bot- 
tom of the band at (tt/2, tt/2, tt/2) agrees with the SDW MF 
result. The QMC data, also shown, lie between these weak 
and strong coupling approximations. The area in the circles 
is proportional to the peak intensity. The error bars are the 
width of the peak. 

Actually, the dispersion of a dressed hole in an anti- 
ferromagnet within the SCBA for a bilayer system, and 
also for a 3D cubic lattice, has been recently studied. pOjl 
Here, for completeness, we reproduce some of the results 



of Ref . |3!| , and compare them against those of the 3D 
Hubbard model obtained with the SDW MF approxima- 
tion and QMC calculations (Fig. 6). The comparison is 
carried out at J/t ~ 0.3 which corresponds to U/t ~ 13. 
The maximum of the dispersion in the valence band us- 
ing the SCBA now lies at (tt/2, tt/2, tt/2), removing the 
spurious SDW MF degeneracy. In the scale of Fig. 6 
the splitting between this momentum and (tt, 7r/2, 0) is 
difficult to resolve, since it corresponds to about 100K. 
Note that the bandwidth predicted by the SDW MF 
technique is approximately a factor two larger than the 
more accurate prediction of the SCBA. However, for this 
larger value of U it does not appear possible to fit si- 
multaneously the SDW MF bandwidth and band-gap 
to the results of QMC by the same renormalization of 
U, something which can be done successfully at weaker 
coupling, U = 8. The QMC points at this intermedi- 
ate coupling value where t/=bandwidth lie in between 
the SDW MF and SCBA. Though the uncertainties in 
the QMC results are rather large, we expect the agree- 
ment between SCBA and QMC results to improve as 
the coupling increases. The best fit of the SCBA data 
p9| is e(p) = c + 0.082(cosp x cosp y + cosp y cosp z + 
cosp x cosp z ) + 0.022(cos2p x + cos2p a + cos2p z ) (eV), 
if J = 0.125eV and t = OAeV are used. The constant 
c is defined by the SDW MF gap (Fig. 6). As in the 
case of the 2D problem, holes tend to move within the 
same sublattice to avoid distorting the AF background. 
Q Working at small J / t, the bandwidth of the 3D t — J 
hole quasiparticle was found to scale as J, ||)) as occurs 
in two dimensions. 



IV. FINITE HOLE DENSITY 
A. D=3 

We can also use the QMC approach to study the 3D 
Hubbard model away from half-filling for temperatures 
down to about 1/30 of the bandwidth, a value for which 
T ~ J for the present strong coupling values. First, 
we study the influence of doping and temperature on 
the spin-spin correlation function C(l). At half filling 
C (1) shows strong antiferromagnetic correlations over the 
whole 4 3 -lattice at ft = 10 (Fig. 7). At = 2 the cor- 
relations are significantly weakened, and with additional 
doping ((n) — 0.88) all correlations are suppressed be- 
sides those between nearest neighbors. These appear to 
be stab le against doping. The density of local moments, 
y C(0) reaches its low temperature limit at an energy 
scale set by U and hence is unaffected by the change of 
(3 from fi = 2 to (3 = 10 (note that longer range spin cor- 
relations form at a temperature set by the much smaller 
energy scale J). y/C(0) is to first order proportional 
to the electronic density and hence slightly reduced at 
(n) = 0.88. 

There has been considerable discussion concerning the 



G 



relationship between the spin-spin correlations and the 
presence of a gap in the density of states. In particular, 
it was observed jll| that if N(u) is evaluated on lattices 
of increasing size at fixed temperature a well formed gap 
appearing on small lattices disappears when the spatial 
extent exceeds the spin-spin correlation length. Decreas- 
ing the temperature (and hence increasing the range of 
the spin correlation) allows the gap to reform. Similar 
effects are seen here in 3D. 



0.3. 
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FIG. 7. Spin-spin correlation function C(l) over a path in 
real space on the 4 3 lattice at U = 8, /? = 10, (n) = 1 (o); 
P = 2, (n) = 1 (x); and j3 = 2, (n) = 0.88 (o). C(0)-values 
are divided by a factor of 3 for clearness. Error bars are 
smaller than the symbols when not shown. 

Fig. 8a shows the density of states on a 4 3 lattice at 
several densities, U = 8 and (3 = 2. At this temperature 
the charge gap is not fully developed, and the quasipar- 
ticle peaks cannot be resolved. The result with doping is 
similar to that reported on 2D lattices. The chemical 
potential /x moves to the top of the valence band as the 
density is reduced from half-filling. A large peak is gener- 
ated which increases in intensity as (n) is further reduced. 
The weight of the upper part of the spectrum (reminis- 
cent of the UHB) decreases with doping due to the re- 
duced effective interaction. Similar results are shown in 
Fig. 8b but for a 6 3 lattice. There is not much difference 
between the two lattices, showing that within the resolu- 
tion of the ME procedure finite size effects are small. 

The large peak that appears in Fig. 8a-b at finite hole 
density is crossed by /j, as the density is reduced. At 
(n) = 0.94, the peak is located to the left of /i, at 
(ri) = 0.88 it has reached the chemical potential, and 
at (n) — 0.72, the peak has moved to the right. This 
is in agreement with the behavior observed in both 2D 
QMC and ED simulations, jllj and it may be of relevance 
for estimations of superconducting critical temperatures 
if a source of carrier attraction is identified. p2| 

The results of the previous section at half-filling ob- 
tained at low temperatures (T ~ 1/10) revealed a sharp 



quasiparticle peak in the DOS at the top of the valence 
band and bottom of the conduction band. Numerical 
studies of 2D lattices have shown that the peak inten- 
sity at T = is the largest at half-filling. |^| Away 
from half-filling, the peak is still visible but it is broader 
than at (n) — 1. [fll] Thus there is no evidence that the 
sharp peak in the DOS of the doped system has been gen- 
erated dynamically and represents a "Kondo-resonance" 
induced by doping, as has sometimes been suggested |L2] , 
and as Fig. 8 obtained at relatively high temperature, 
(3 = 2, seem to imply. 
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FIG. 8. (a) N(lo) for the 3D Hubbard model on a 4 3 lattice 
with (7 = 8 and ft = 2 at several densities indicated in the 
figure; (b) same as (a) but on a 6 3 lattice. Frequencies are 
always relative to the chemical potential [i. 

Another important quantity to study is the quasipar- 
ticle residue Z. The SCBA results show that Z is small 
but finite for the case of one hole in an antiferromagnetic 
insulator state, and actually the results are very similar 
in 3D and 2D systems. |38|;p9| Numerical results provide 
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a similar picture, g] On the other hand, Z vanishes in the 
D = oo approach working in the paramagnetic state as 
the doping 5 tends to zero. Note that in this state there 
are no AF correlations (£af = 0). Thus, it is clear that 
the hole quasiparticle at half-filling observed in the 2D 
and 3D systems is not related with the quasiparticle-like 
feature observed in the PM state at D = oo. 
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FIG. 9. (a) Spectral weight A(p,co) of the 3D Hubbard 
model calculated using QMC/ME on a 4 3 lattice, at U = 8, 
j3 = 2, and at density (n) = 0.94; (b) same as (a) but at 
(n) = 0.88; (c) same as (a) but at (n) = 0.72. 

In Fig. 9, we show A(p, u>) obtained on the 4 3 lattice, 
U = 8, T = 1/2 and various densities away from half- 
filling. The gap is now absent. From the energy location 
of the maximum of the dominant peak in Fig. 9a-c, the 
quasiparticle dispersion can be obtained. The results are 



shown in Fig. 10. It is remarkable that the quasiparticle 
dispersion resembles that of a noninteracting system i.e. 
e p = —2t*(cosp x + cos p y + cosp z ), with a scale increas- 
ing from t* ~ t/4 to i/3 with doping. This dispersion 
certainly does not exhaust all the spectral weigth but a 
large incoherent part still remains at this coupling, den- 
sity and temperature. Similar results were observed in 
2D. [p] j2^ , [I]] . [l3|] Only vestiges remain of the AF induced 
weight in PES near (tt, 7T,7t). However, this drastic re- 
duction of the AF induced intensity may be caused by 
the high temperature of the simulation as observed in 
the spin-spin correlation function (Fig. 7). 
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FIG. 10. Dispersion of the dominant peak of Fig. 9a-c 
against momentum. The densities are indicated. The area in 
the circles is proportional to the peak intensity. Error-bars 
correspond to the half-width of the peaks in the spectral 
weight. 



B. D = oo 

The previous subsection and the results at half-filling 
have shown that the DOS of the 3D Hubbard model has 
a large peak at the top of the valence band. The peak 
is crossed by the chemical potential as (n) decreases. 
This behavior is in apparent contradiction with results re- 
ported at D = oo where a peak is generated upon doping 
if the "paramagnetic" solution to the mean-field problem 
is selected. At D — oo, there are only two very distinct 
magnetic ground states. One has AF long-range order, 
and the other is a paramagnet with strictly zero AF cor- 
relation length i.e. without short range antiferromagnetic 
fluctuations. Thus, at D — oo the transition is abrupt 
from a regime with £af = oo to £af = 0. This does not 
occur in finite dimensions where before the long-range 
order regime is reached, AF correlations start building 
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up smoothly. This qualitative difference is depicted in 
Fig. 11. 
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FIG. 11. "Phase diagram" of the Hubbard model (1) in 
the D-(n)— plane. Solid line: AF phase boundary; dashed 
line: crossover where short range AF correlations disappear. 
The intermediate regime of short range AF order vanishes in 
the limit D — > 00. 

S;af as small as a couple of lattice spacings can be ro- 
bust enough to induce important changes in the carrier 
dispersion, and may even be enough to induce supercon- 
ductivity as many theories for the 2D high-Tc cuprates 
conjecture. We believe that the absence of a regime of 
intermediate size AF correlations at large D is the key 
ingredient that explains the differences reported here be- 
tween D=2,3 and D = 00. 

In Fig. 12a, the D=oo DOS in the AF phase is shown 
at (n) = 1 and 0.94. For these densities the AF-phase 
is energetically stable. We observe the tendency of the 
large peak at the bottom of the valence band to move to- 
wards the chemical potential in good agreement with the 
3D Quantum Monte Carlo simulations. As found in 2D, 
the intensity of the peak decreases as we move away from 
half-filling if the temperature is low enough. In Fig. 12b, 
the DOS in the D = 00 limit working in the paramag- 
netic phase is shown at several densities. For the present 
interaction, U = 8, the paramagnetic solution remains 
metallic at all temperatures even at half filling. The 
results are qualitatively different from those observed in 
the AF regime. At (n) = 1 a large peak at the chemical 
potential is clearly visible. Upon hole doping this peak 
gradually moves toward higher energies. At sufficiently 
strong doping the DOS of the PM-phase (Fig. 12b) re- 
sembles the results for the 3D lattices (Fig. 8), which is 
not surprising since AF correlations in 3D are strongly 
suppressed at the present temperature. Close to half fill- 
ing, however, the 3D results are closer to the DOS of 
the AF-phase where a strong peak is observed on the 
left hand side of the chemical potential /i. This result is 



gratifying since the proper way to compare D = 3 and 
00 results is by using the actual ground states in each 
dimension. In D = 00, at low temperatures, the cross- 
ing of the peak by fj, is expected at that point where the 
AF-Phase becomes unstable against doping. 



0.20 - a) 



0.15 - 



n=1.00 
n=0.94 



3 

Zf 0.10 



0.05 - 




U=8 
P=4 
AF 



0.00 1 1 1 

-12-8 -4 4 

CO 



0.18 




8 12 



s 



0.12 - 



0.06 



0.00 



b) 




n=1.00 


U=8 


f 93 


\ p=2 


n=0.87 

jl 
If 


\ PM 




\ i \ v 



-15 



-10 





CO 



10 



15 



FIG. 12. (a) N(uj) corresponding to the D = 00 Hubbard 
model at U = 8, /3 = 4 at the electronic densities indi- 
cated. The results were obtained using the AF solution to 
the mean-field equation; (b) same as (a) but using (3 = 2 and 
the PM solution to the mean-field equation. 



V. CONCLUSIONS 

In this paper we have calculated the single particle 
properties of the 3D single band Hubbard model using 
Quantum Monte Carlo and the SDW mean field and 
SCBA approximations. Our results have many similar- 
ities with those reported previously in 2D systems. At 
half-filling, peaks at the top of the valence band and bot- 
tom of the conduction band are observed in the DOS. 
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Their behavior is associated with spin polarons with a 
bandwidth of order the exchange J. We found similar- 
ities to and semi-quantitive agreement with experimen- 
tally observed features in the spectra of strongly corre- 
lated 3D AF insulators, LaFeC>3 and V2O3. 

As we dope the system, the sharp peak associated with 
these quasiparticles is crossed by the chemical potential 
as the density (n) changes. The PES weight observed 
away from half-filling is already present at half-filling. 
No new states are generated by doping. This result must 
be contrasted with that observed experimentally in, e.g., 
Yi_ x Ca x Ti03 using angle-integrated PES. In this case 
spectral weight which is not present in the insulator ap- 
pears at Ep in the metallic regime as we dope the sys- 
tem. This behavior does not seem reproduced by the 
single band Hubbard model Eq.(l) in 3D, whose physics 
appears to be very close to that of 2D. Indeed for the 
2D cuprates it has been shown experimentally that the 
states found at Ep upon doping are already present at 
half-filling. @ 

An exception among 3D materials is NiS2_ x Se x 
which remains antiferromagnetic throughout the metal- 
insulator transition induced by (homovalent) Se- 
substitution or temperature. PES spectra at x = 0.5 
for different temperatures show a strong peak close 
to the Fermi energy which does not disappear in the in- 
sulator. Instead the peak is shifted off the Fermi energy 
and only very slightly reduced in weight. Since this sit- 
uation is not described within the paramagnetic D = 00 
approach, AF correlations are presumably essential for 
the low energy electronic excitations of this system. 

The success of the D = 00 approach to the Hubbard 
model in describing the physics of Yi_ x Ca x Ti03, SrVC>3 
and CaVC>3, however, appears crucially to depend upon 
forcing the paramagnetic solution of the equations. |4(| 
In this case, states are actually generated in the Hubbard 
gap after a small hole doping is introduced. Of course, it 
may be that the "arbitrary" choice of this paramagnetic 
solution, which is not the actual minimum of the free 
energy, is well motivated since it mimics the presence of 
physical effects like frustration which destroy long range 
order in real materials. More work is needed to show 
that this scenario is realized for realistic densities and 
couplings. 

An alternative explanation for the discrepancy between 
the PM solution in infinite-D and finite dimensional re- 
sults may lie in the finite resolution of the combination 
of Monte Carlo simulations and Maximum Entropy tech- 
niques. However, the SCBA and results at half-filling 
and low T show that it is likely that at (n) = 1 we have 
quasiparticle states in the DOS. 

In studies of the single band Hubbard Hamiltonian in 
2D, and in the present analysis in 3D, it is clear that 
short-range AF correlations play an important role close 
to (n) = 1. In particular, the states created at the top 
of the valence band are likely to be spin polarons with 
a finite quasiparticle residue Z. PES states observed at 
finite hole doping evolve continuously from those present 



at half-filling. Experiments on the 2D high-Tc cuprates 
seem to present similar features, while the results for the 
3D perovskites are very different in the sense that no 
remnants of the coherent part of the spectrum away from 
half-filling are reported at half-filling. 

Still, strong AF correlations are apparently present in 
several 3D transition metal oxides and influence the low 
energy spectrum at least on the insulating side of the 
transition. The introduction of frustration in the sin- 
gle band Hubbard model in 3D, perhaps through next- 
nearest neighbor hoppings, will reduce AF correlations 
and in particular the AF-induced charge gap, and might 
be sufficient to observe a evolution of spectral weight 
upon doping closer to the experimental findings. How- 
ever, it might be that 3D models which explicitly in- 
clude orbital degeneracy will be necessary to reproduce 
the physics of the transition metal oxides, as has recently 
been described for NiO chains |4^| and Mn-oxides. [ EgJ . 
Indeed, a recent argument presented by Kajueter et al. 
p9[ to justify the use of the D — 00 model provides a 
more realistic explanation for the apparent link between 
theory in this limit, and 3D transition-metal oxide re- 
sults. The idea is that the physics of the real perovskite 
3D oxides is influenced by the orbital degeneracy. Pre- 
sumably this effect leads to a drastic reduction of the 
antiferromagnetic correlations that dominate the physics 
of these 2D and 3D systems. Many orbitals, including 
Hund's coupling, produce an effective magnetic frustra- 
tion that may reduce the AF correlation length to a neg- 
ligible value even close to the AF insulator at half-filling. 
Such a frustration effect could be strong enough to gen- 
erate a finite critical coupling U/t at half-filling. 
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